Quantum Phase Transitions and Continuous Observation of Spinor Dynamics in an 

Antiferromagnetic Condensate 
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Condensates of spin-1 sodium display rich spin dynamics due to the antiferromagnetic nature of 
the interactions in this system. We use Faraday rotation spectroscopy to make a continuous and 
minimally destructive measurement of the dynamics over multiple spin oscillations on a single evolv- 
ing condensate. This method provides a sharp signature to locate a magnetically tuned separatrix 
in phase space which depends on the net magnetization. We also observe a phase transition from a 
two- to a three-component condensate at a low but finite temperature using a Stern-Gerlach imaging 
technique. This transition should be preserved as a zero-temperature quantum phase transition. 

PACS numbers: 67.85.-d, 03.75.Kk, 03.75.Mn, 64.70. Tg 



The study of multi-component Bose-condensed (super- 
fluid) systems began with ^He - ®He mixtures in the 
1950's |l| and continued with a two-component *^Rb 
Bose-Einstein condensate (BEC), with atoms in two dif- 
ferent hyperfine states, in the late 1990's 2] . Recent stud- 
ies have focused on the investigation of spinor conden- 
sates where interconversion among multiple spin states 
leads to spin population dynamics. A number of investi- 
gations of this effect, as well as spatial domain formation 
in both ferromagnetic F=l *^'^Rb BECs 0,0, H and an- 
tiferromagnetic F=l 2^Na BECs [i,|3,[l have been pub- 
lished. The F=2 s^Rb spinor condensate presents ferro- 
magnetic, antiferromagnetic, and cyclic phases 
H^ . In each of these cases the experimental system can 
be modelled with a small number of variables. 

While the ferromagnetic Rb system is becoming well- 
studied, the antiferromagnetic Na system remains rela- 
tively unexplored. The theoretical description of such a 
system suggests the possibility of manipulating the phase 
space topology and dynamics of the system in ways not 
possible in the ferromagnetic system (for example, al- 
tering the separatrix position in phase space with the 
magnetization of the system) as well as the possibility of 
observing a quantum-fluctuation-driven phase transition 
that does not exist in the ferromagnetic system [1, . 

We consider a spinor BEC with spin angular momen- 
tum = 1 in the presence of a magnetic field of strength 
B along the z axis with the populations initialized to 
a non-equilibrium state. CoUisional interconversion be- 
tween two mp = Q atoms and one mp — +^ and one 
nip = — 1 atom takes place in the condensate, leading to 
oscillations in the spin populations. At ultracold temper- 
atures the collisions between alkali metal atoms conserve 
the summed spin angular momentum f = Fa + Fi,. Our 
system, ^^Na, is antiferromagnetic, or polar, inasmuch as 
the interaction energy of / = 2 collisions is larger than 
that of / = collisions, which indicates that the coupling 
favors the mp = ±1 states over the mp — state. 



The linear Zeeman shift induced by the magnetic field 
does not affect the collisional interconversion, as the mag- 
netic energies before and after the collision are equal in 
this approximation. The population dynamics are in- 
stead driven by an interplay between the quadratic Zee- 
man shift and the spin-dependent interaction character- 
ized by the difference in the / = and / = 2 interaction 
energies. In ^^Na, a divergence in the spin oscillation 
period occurs near a critical magnetic field Be @, [sl, [l^ . 
A dependence of Be on magnetization m (the difference 
in fractional population p^p between mp — +1 and 
mp = —1) is predicted in an antiferromagnetic system 
or in a ferromagnetic system with a radio-frequency (rf) 
dressing field [14|, but has not been previously observed. 

We use two complementary methods to observe the 
spin dynamics in two different time regimes. A Stern- 
Gerlach separation followed by absorption imaging (SG- 
AI) is the standard method to directly measure the pop- 
ulations of different spin states and determine the mag- 
netization m. This technique is, however, completely de- 
structive and only minimal phase information about the 
spin oscillations can be inferred from modelling the data. 
The second method is Faraday rotation spectroscopy 
which measures the rotation of the polarization of a laser 
beam. This rotation is proportional to the projection 
of the atomic spin F along the laser propagation direc- 
tion. It can be used to continuously infer both phase and 
population information of the spin dynamics over mul- 
tiple spin oscillations. Other methods of measuring the 
condensate phase can be found in 
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It is hard to determine B^ from just the spin oscillation 
period, however, we observe a sharp signature to distin- 
guish two characteristic time evolutions in the vicinity of 
Be with Faraday rotation spectroscopy. At long times 
when the oscillations have damped out 0], we use SG-AI 
to characterize the mean-field ground state populations 
as a function of B and m to observe a predicted phase 
transition from a two- to a three-component spinor BEC. 
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A ^^Na magneto-optical trap containing up to 6 x 10^ 
atoms is prepared. A crossed optical dipole trap de- 
rived from a multi-mode laser at 1070 nm is then loaded, 
followed by evaporation and re-thermalization. A weak 
magnetic field gradient is applied during 6 s forced evap- 
oration to form a fully polarized BEC of 1.5 x 10^ atoms 
in the \F = l^mp = 1) state. The final trap frequencies 
are lOx y z « (27r)220(V^, 1, 1) Hz, and the mean Thomas- 
Fermi radius is 7.2 /im. We ramp up the magnetic field 
along the z axis while turning off the field gradient. The 
final value of B ranges between 6.3 /xT and 60.7 /iT 
with an uncertainty of 0.04 /iT (uncertainties in this Let- 
ter represent estimated one standard deviation combined 
statistical and systematic uncertainties) . We can prepare 
an initial state with any desired m and pq by an rf-pulse 
to rotate the atomic spin followed by selective removal 
of atoms in a given mp state. The rf-pulse is resonant 
with the linear Zeeman splitting, and its amplitude and 
duration control the superposition of the mp levels. The 
removal is performed by a microwave pulse to selectively 
transfer \F = l^mp) atoms to the F — 2 state, followed 
by a laser pulse resonant with these atoms. 

The Faraday detection beam is directed along the x 
axis and red-detuned 225 GHz from the 10 MHz wide 
D2 line of ^'^Na. The beam is linearly polarized, has a 
1/e^-waist of 1 mm at the condensate, and a power of 
~50 mW. The set-up for Faraday spectroscopy is similar 
to that outlined in 1^ 17 1. A carefully chosen aperture 



is inserted into the imaging plane for an optimal signal- 
to-noise ratio (SNR), and the Faraday rotation of the 
linear polarization is detected using a WoUaston prism 
and an autobalanced differential photodetector (PD) [3]. 
The Faraday rotation angle oscillates at the Larmor pre- 
cession frequency Jl- Changes of spin populations and 
phases are detected as a modulation of the amplitude of 
the Faraday signal. Our Faraday signal is the short-time 
power spectral density of the PD output integrated over 
a narrow bandwidth of 1 kHz around /l . This is propor- 
tional to the slowly varying envelope of {F^)'^ . A typical 
example is shown in Fig.[TJ We divide the Faraday signal 
into 1 ms time bins, longer than the transform limit of 
the digital filter (0.16 ms) but short enough to resolve 
spin oscillations. Over a 100 ms measurement, we thus 
make 100 distinct measurements of both the spin projec- 
tion amplitude and /i, on a single evolving spinor BEC. 

The parameters of the Faraday beam are chosen to 
minimize atom loss from the off-resonant light scattering 
while maintaining a good single-shot signal. The mea- 
sured lifetime of the BEC in the presence of the Fara- 
day beam is 100 ms, consistent with the decay of {F^) 
inferred from Fig. [U and with the predicted photon scat- 
terin g ti me. The dephasing time due to the tensor light 
shift [19| is one order of magnitude longer. The scattering 
loss is larger than any other back-action in the present 
experiments. In the absence of the Faraday beam, we 
observe an energy dissipation which depends on B and 



the mean particle density (n). Under the conditions of 
Fig. [ija), the time scale of this dissipation is five times 
longer than the decay seen in this figure, while at high 
fields (Fig. [Ijb)), it becomes comparable. This dissipa- 
tion is not well understood. 

The SNR of our measurements is limited by the total 
number of atoms in the BEC and the efficiency of the 
detection [3] ■ Our BECs are not much larger than what 
is required to get a good single-shot signal with our sys- 
tem. Our overall detection efficiency could be improved 
by a factor of two. 
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FIG. 1; Faraday signal (proportional to (Fx)^) taken from a 
single measurement for m = at two magnetic fields, 26 fiT 
and 40 /iT starting with po ~ 0.5, 6 = 0. The solid line is a 
fit with a damped sinusoid. The signals show an oscillating 
phase (a) and a running phase (b) at B below and above Be, 
respectively, as evidenced by the signal reaching zero or not. 

The single mode approximation (SMA) [^, @] is applied 
to understand our data. The spatial wavefunction of the 
BEC is treated as a single mode and the unit-normalized 
total wavefunction can be represented as '^{r,t) = 
<I'(r)(V^r7(Z)e^«-(*), y;^e^«»(*), VmWe^'+^^*^),^ 
where and phases are independent of position. 
The Hamiltonian conserves particle number and m. The 
system is described using po and 6 = 0-i + 0+i — 29q, 
with the conserved classical spinor energy 

E = £;qz(l-po)+cpo ((1 - po) + pcf -m'^ cos0 

(1) 

Here Eqz is the quadratic Zeeman shift (Eq^/h — 
(0.0277Hz/(pT)2)B2), c = C2{n) is the spin-dependent 
collision energy, and C2 is 1.59 x 10~^^ Jm'^ for ^^Na [1]. 
The evolution of po and 6 is given by po = —{2/H)dE/d9 
and e = {2/h)dE/dpQ. 

Figure O shows typical phase diagrams of the equal 
energy contours of Eq. [T] for m = at two magnetic 
fields. The preparation of the state determines the energy 
Eo{B). Our initial states have 6 = 0. At any magnetic 
field, we can define a separatrix, i.e. that energy contour, 
EscpiB), on which there is a saddle point where po = 
6 = 0. This defines the boundary between two regions in 
phase space. In fact, Escp{B) = Eqz in our system. 
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When Eq{B) > Escp{B), the value of 9 is restricted, 
while for Eq{B) < Escp{B) there is no bound. This de- 
fines regions with oscillating phase and running phase, 
respectively. In both regions po oscillates. At Eq{B) ^ 
Esep, which corresponds to small magnetic fields, the pe- 
riod only weakly depends on the field. In the opposite 
limit, the period decreases rapidly with increasing field. 
In both cases, the oscillation is nearly harmonic. When 
B ^ Be, defined by Eq{Bc) = Escp{Bc), the oscillation is 
anharmonic and the period diverges for B = Be- 



B = 26 |jT 
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FIG. 2: (Color online) Equal-energy contour plots generated 
from Eq. [1] at two magnetic fields, 26 (left) and 40 /^T 
(right), with m = and c/h — 33 Hz. The (red) dashed lines 
represent the energy of a state with po{t = 0) = 0.5. The 
heavy (blue) solid lines represent the energy of the separa- 
trix (-Escp) between oscillating and running phase solutions. 
Darker colors represent lower energies. 

In the SMA, the Faraday signal is derived from 
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y/po (1 + m - Po) (2) 
Vpo (f - m- Po). 



The phase difference 6+i — O-i is determined by the fast 
Larmor precession and a slow evolution due to po and 
9 For m = 0, our Faraday signal is proportional to 
/Oo(f — Po) cos^(6'/2). For oscillating phase solutions, 9 
oscillates about zero (with amplitude < tt) and thus the 
signal is always greater than zero. On the other hand, 
the signal is periodically zero for running phase solutions. 
Figure [T] shows signals from the two regions. For m ^ 0, 
the signal is described by a more complicated expression, 
but the distinction between the two regions in the vicinity 
of Be remains the same. 

These characteristics of the Faraday signal provide a 
sharp signature for locating Be- Figure [3] shows the value 
of the minimum of the Faraday signal at different B and 
771 after removal of the exponential decay. For the two 
magnetizations a transition from an oscillating phase so- 
lution with a non-zero minimum to a running phase so- 
lution with a minimum of zero, provides the sharp signa- 
ture to locate Be- In the population oscillations were 



measured using SG-AI and fit by a sinusoid to extract a 
period and to locate the two regions of phase space. No 
sharp experimental signature distinguishing the bound- 
ary was observed. 
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FIG. 3: (Color online) The minimum of the Faraday signal 
as a function of magnetic field for m = (red dots) and m = 
0.3 (blue bowties). A scale factor is applied to the Faraday 
signal to correct for the PD response at different /l. The 
lines are fits based on Eq. [2] which yield the fit parameters 
Po = 0.42 and iV = 1.50 x 10^ for m = 0, and po = 0.54 and 
= 1.32 X lO'' for m = 0.3. The fit parameters are within 
the 3% uncertainty of those derived from absorption images. 

Figure [3] also shows a comparison between the pre- 
diction from the SMA and the data. For 771 = 0, the 
agreement is excellent. For 777 — 0.3, however, the pre- 
diction does not agree with our measurements for fields 
significantly larger than B^- At the transition point, 
the minimum of the Faraday signal goes to zero, but 
above this point, the theory predicts that the minimum 
rises with B, even though the solution has a running 
phase. This increase is not observed. The apparent 
agreement between the SMA theory and our measure- 
ments has been surprising - at every magnetic field re- 
ported in this paper, we have seen the presence of several 
spatial modes/spin domains during the spin oscillations, 
although not in steady-state. The observation of spin 
domains is in marked contrast to our previous work [8(]. 
Several technical changes may have contributed to the 
domain formation, such as a 50% increase of the total 
atom number and more stable magnetic bias fields. Un- 
derstanding the ways in which multiple spatial modes 
and domains affect the spin dynamics is an interesting 
direction for future research. 

At long times, after the oscillations have damped out, 
we can study the mean-field ground state of this sys- 
tem. Within the SMA, a quantum phase transition from 
a three- to a two-component BEG is predicted for the 
mean- field ground state [l3[. We use SG-AI as a direct 
way to measure the equilibrium populations and to ob- 
serve this phase transition. 
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B = 21.2 [iT 
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FIG. 4: Evidence of phase transitions in the mean-field 
ground state of the antiferromagnetic spinor BEG. (a): Points 
with uncertainties indicate po as a function of _B at m = 
0.71(2). The intersection of two linear fits defines -B2-3, the 
critical magnetic field for the phase transition, (b) Inset: po 
versus m for B — 21.2 fj,T. The intersection of two linear 
fits defines 7712-3, the critical magnetization. The main figure 
shows m2-3 versus B. The solid line is the prediction from 
the SMA. 



prediction from the SMA as shown in Fig. [H^b). 

This confirms a phase transition from a three- to a two- 
component spinor BEC in the antiferromagnetic mean- 
field ground state. Although it is observed here at a finite 
temperature the phase transition should remain at zero 
temperature, where the transition would be driven solely 
by quantum fluctuations. It has been predicted that this 
phase transition persists even when the three spin states 
do not share the same spatial distribution and the appli- 
cation of the SMA is no longer appropriate • 

In conclusion, we have demonstrated that Faraday ro- 
tation spectroscopy provides a method to continuously 
monitor the spin dynamics in an antiferromagnetic spinor 
BEC. The technique provides a sharp signature to locate 
the magnetically tuned boundary in phase space between 
the oscillating and running phase solutions. We observe 
a dependence of the position of this separatrix on the 
magnetization. In addition, we have confirmed a quan- 
tum phase transition from a two- to a three-component 
BEC in the mean-field ground state. We are presently 
investigating the dissipation mechanism leading to the 
equilibrium observed in these experiments. We think 
that physics beyond the SMA, possibly similar to Landau 
damping of excitations of a single component BEC 2l|, 
may explain the dissipation. 

We thank W. D. Phillips, J. V. Porto and E. Gomez 
for insightful discussions, and ONR for financial support. 
SEM thanks the NIST/NRC postdoctoral program. 



The behavior of the variance of po allows us to deter- 
mine when the system has settled down to a steady state. 
At a given time we measure po 25-30 times and calculate 
a variance. In the steady state, the variance reaches a 
minimum. Due to technical noise in atom counting, this 
variance is larger than that predicted by a quantum cal- 
culation of the spinor ground state [20]. We find that a 
steady state for to = is reached within our maximum 
hold time of 10 s for B > 18 p,T. For non-zero m, the 
steady state is reached within a much shorter time. 

Figure S] (a) shows the steady state values of po as a 
function of B. For each measurement, we prepare nearly 
identical initial states a,t B = 24.4 ^T to set the popu- 
lations and magnetizations. We then wait 4 s to reach 
the steady state before adiabatically ramping the field to 
a desired final value over 500 ms. We then wait another 
5s before making an SG-AI measurement. In the inset of 
Fig. 131(b), however, each initial state is directly prepared 
with a different magnetization. In both figures, a tran- 
sition between a two-component BEC with pQ ~ and 
a three-component BEC with non-zero po, is observed. 
A critical magnetic field B2-3 and a critical magnetiza- 
tion 7712-3 are defined and extracted from the intersec- 
tion of two linear fits to the data. Good agreement is 
found between the experimental value of TO2-3 and the 
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